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ABSTRACT 

We draw attention to the close relationship between the anisotropy parameter (3{r) for stellar orbits in 
elliptical galaxies and the temperature T(r) of the hot interstellar gas. For nearly spherical galaxies, the 
gas density p can be accurately determined from X-ray observations and the stellar luminosity density 

can be accurately found from the optical surface brightness. The Jeans equation and hydrostatic 
equilibrium establish a connection between [3{r) and T(r) that must be consistent with the observed 
stellar velocity dispersion. Optical observations of the bright elliptical galaxy NGC 4472 indicate /3(r) < 
0.35 within the effective radius. However, the X-ray gas temperature profile T(r) for NGC 4472 requires 
significantly larger anisotropy, /3 « 0.6 — 0.7, about twice as large as the optical value. This strong 
preference for radial stellar orbits must be understood in terms of the formation history of massive 
elliptical galaxies. Conversely, if the smaller, optically determined (3{r) is indeed correct, we are led to 
the important conclusion that the temperature profile T{r) of the hot interstellar gas in NGC 4472 must 
differ from that indicated by X-ray observations, or that the hot gas is not in hydrostatic equilibrium. 

Subject headings: galaxies: elliptical and lenticular, CD - galaxies: active - cooling fiows - X-rays: 
galaxies - galaxies: clusters: general - X-rays: galaxies: clusters 



1. introduction 

The anisotropy of stellar orbits in elliptical galaxies can 
be estimated from the temperature of the hot interstel- 
lar gas through which they move. These two quite differ- 
ent galactic attributes are intimately related by the Jeans 
equation for the stars and the condition for hydrostatic 
equilibrium in the gas. This simple relationship is valu- 
able since both the anisotropy and gas temperature are 
difficult to extract from optical and X-ray observations, 
respectively. 

Many massive elliptical galaxies are nearly spherical 
(Merritt & Trembly 1996) and slowly rotating, but their 
stellar velocity ellipsoids are not in general isotropic. The 
non-spherical nature of the stellar velocity dispersion is 
represented by the parameter 

/3 = 1 - aVal 

where Ur is the radial stellar velocity dispersion and at is 
the dispersion in a transverse direction, i.e. — — a"^. 
If the orbits are predominantly radial, < /3 < 1, the line 
of sight velocity profile becomes more strongly peaked than 
a Gaussian profile with increasing projected radius R\ if 
the orbits are mostly tangential, — oo < < 0, the profile 
is more flat-topped and becomes broader with increasing 
radius. Both /3(r) and the galactic potential $(r) can be 
determined from optical observations of the line of sight 
velocity dispersion cr(-R), the optical surface brightness dis- 
tribution and the deviation of the stellar line profiles from 
a Gaussian as expressed by the line-symmetric coefficient 
hi{R) in a Gauss-Hermite expansion (e,g, van der Marel 
& Franx 1993). 

Stellar line profiles observed in most massive E galaxies 



indicate a preference for radial orbits with j3 ~ 0.3 (Ben- 
der, Sagha & Gerhard 1994; Gerhard et al. 1998; Sagha et 
al. 2000). Accurate /3(r) profiles from optical data require 
high quality data and considerable care in reduction and 
analysis. The radially anisotropic nature of stellar orbits in 
luminous elliptical galaxies provides important and other- 
wise unavailable information about the merger history of 
these galaxies (Naab, Burkert & Hcrnquist 1999), so an 
improved or independent determination of [3[r) would be 
desirable. 

Similarly, the radial variation of the hot gas tempera- 
ture T{r) ~ Tyir ~ 10^ K in elliptical galaxies depends on 
the spatial and spectral resolutions of X-ray detectors and 
the accuracy of the three-dimensional decomposition that 
converts T as a function of projected radius R to physical 
radius r. In addition, it is often unclear whether the hot 
gas at any radius has a single temperature or a multitude 
of temperatures as might be expected if the gas were cool- 
ing. It is unclear if the gas cools at all. For example, high 
resolution X-ray spectra with XMM-Newton of the large 
E galaxy NGC 4636 fail to show emission lines expected 
from gas at intermediate temperatures, such as the 0.574 
keV OVII line (Xu et al. 2002), suggesting that the gas 
is not cooling below ~ Tyir/3 as in classical cooling flows. 
However, Bregman, Miller & Irwin (2001) detected the 
OVI 1032,1038A doublet in NGC 4636 emitted from gas 
at r ^ 3 X 10^ K, implying that cooling near the expected 
rate may occur after all. As an additional source of confu- 
sion, the observed X-ray spectra can often be signiflcantly 
improved by assuming two quite different discrete temper- 
atures at each galactic radius (e.g. Buote 2002 for NGC 
1399; Buote et al. 2002a and Tamura et al. 2003 for NGC 
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5044). While the origin and physical nature of gas with 
only two (or a limited continuum of) temperature phases 
may be conceptually problematical, it may be related to 
certain heating mechanisms. Nevertheless, there are good 
reasons for suspecting the acciiracy or interpretation of gas 
temperature profiles T{r) in elliptical galaxies. 

In the following discussion we illustrate the close re- 
lationship between the radial variations of /3(r) and T(r) 
using data for the well-observed giant elliptical NGC 4472. 
While currently available data for this galaxy are adequate 
to illustrate this relationship, more accurate optical and 
X-ray data will be necessary to verify /3(r) or T(r) with 
confidence. 

We recognize that NGC 4472 is not perfectly spherical; 
its optical image has been variously classified as El or E2, 
However, optical determinations of /?(r) for NGC 4472 and 
other nearly-spherical elliptical galaxies have been made 
under the assumption of spherical symmetry (e.g. Kro- 
nawitter et al. 2000). Indeed the very concept of the f3 
anisotropy parameter is only valid in spherical geometry. 
Similarly, the gas temperature T(r) determined from X- 
ray observations of NGC 4472 and other nearly spherical 
elliptical galaxies are based on the assumption of spheri- 
cal symmetry. One justification for this is that the stel- 
lar gravitational potential that confines the stars and hot 
gas within the half-light radius is always more spherical 
than the underlying stellar luminosity distribution (e.g. 
Brighenti & Mathews 1996). 

Our discussion of the close relationship between stellar 
anisotropy and gas temperature is a thinly-veiled request 
for observers to concentrate on those luminous E galaxies 
for which the most accurate X-ray and optical data can be 
acquired. 

2. RELATING /3(r) AND t{r) 

The Jeans equation for the radial stellar velocity disper- 
sion ar is 

dr r ^ ^ ' 

where ^* is the stellar luminosity density corresponding to 
the stellar velocity dispersions ar and at- The equation 
for hydrostatic equilibrium in the hot interstellar gas is 

1 d{pc^) 

p dr 

Here c? = kT/fxrUp is the isothermal sound speed, nip is 
the proton mass and jj, = 0.61 is the molecular weight. 
The uniformity of the stellar mass to light ratio within the 

effective radius of NGC 4472 is excellent evidence that the 
hot gas is in hydrostatic equilibrium (Brighenti & Math- 
ews 1997). Both equations contain the same gravitational 

term g = GM(r)/r'^. Eliminating this term results in an 
equation explicitly independent of the confining mass: 



= -9- 



(2) 



da^ 
dr 



dlogi 



d\ogr 



■2f3 



dr 



dlogp 



dlogi 



(3) 



This equation has been discussed previously, particularly 
in the context of the anisotropic orbits of galaxies in rich 
clusters that contain hot gas (e.g. Fabricant et al. 1989). 
In addition to hydrostatic equilibrium, the validity of this 
equation depends on the assumption that the gas pres- 
sure exceeds that of any other interstellar pressure such as 
turbulence, cosmic rays, etc. 



If (3{r) and T(r) (as well as p{r) and i*{r)) are known 
from observations, then Equation (3) can be solved for 
ar{r) and the line of sight velocity dispersion is then found 
from 



a2(i?) 



a^jr) [1 - iR/r)^f^] 4(r)(r2 - R^)-y\dr 
Q i.{r){r^ - R?)-^/^rdr 

(4) 

(e.g. Binney & Mamon 1982). In this equation R is the 
projected radius and rj is some radius significantly be- 
yond the limit of the available data (we assume rt = 3i?e). 
If /3(r) is known. Equation (3) can be solved for (Jr{r) 
for various parameterized temperature profiles T(r) until 
(t{R) from Equation (4) agrees with observed velocity dis- 
persions. By this means it is possible to derive the gas 
temperature profile T(r) from a knowledge of /3(r). Con- 
versely, if the gas temperature T(r) is securely known from 
X-ray observations, then /3(r) can be determined by re- 
quiring that solutions cr^(R) of Equations (3) and (4) pass 
through the observations. This close relationship between 
/3(r) and T(r) is possible and valuable because t*{r), p{r) 
and (T^ (R) are much easier to determine from observations 
than either /3(r) or Tir). 

Except for the very central region, the stellar brightness 
in NGC 4472 is well-fit with a de Vaucouleurs profile. The 
stellar B-band luminosity density ^*(r) = ^*,dev('') is de- 
termined by the total luminosity Lb = 7.89 x 10"'^'' Lb.q 
and the effective radius: Re = 104" = 1.733' = 8.57 kpc, 
assuming a distance d = 17 Mpc. However, within a small 
"break" radius, Vb = 2.41" = 200 pc, the stellar luminos- 
ity flattens to i*^core{T) = (■*,dev{rh){r /rb)~'^ '^^ (Gebhardt 
et al. 1996; Faber et al. 1997). In Figure 1 we com- 
pare the (un- normalized) stellar luminosity density f*(r) 
in NGC 4472 with the hot gas density profile p{r) deter- 
mined from X-ray observations. The remarkable propor- 
tionally ^*{r) oc nf,{r) was first discovered by Trinchieri, 
Fabbiano & Canizares (1986). 

X-ray observations of the hot gas temperature in the 
central region of NGC 4472 are shown in Figure 2. Only 
three observations of the gas temperature are available 
in the region within 10 kpc where (3 has been deter- 
mined from optical observations. The filled circles show 
ROSAT deprojected gas temperature observations from 
Buote (2000). The linear temperature profile defined by 
these two temperature observations 



T(r) = 0.7480 x 10^ + 7.251 x IQ^rkpc 



(5) 



is shown in Figure 2 with a solid line. From Chandra 
data Soldatenkov, Vikhlinin & Pavlinsky (2003) have re- 
cently provided an additional temperature measurement 
of T/10^ K = 0.766 ± 0.023 at radius 1" = 0.0824 kpc 
that is shown with a filled square in Figure 2. The linear 
fit to the Buote data passes within the error bars of the 
Soldatenkov et al. observation, so all three observations 
are consistent with a linear temperature profile within the 
effective radius i?e = 8.57 kpc. Although the X-ray data is 
sparse in the region of interest, it is in perfect qualitative 
agreement with linear gas temperature profiles observed 
in other bright elliptical galaxies in r ^ i?e (Fig- 1 of 
Brighenti & Mathews 1997). 
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2.1. Determining j3{r) from T(r) 

The open circles in Figure 3 show observations of the 
line of sight stellar velocity dispersion in NGC 4472 (from 
Bender et al. 1994) that were used by Kronawitter et al. 
(2000) to determine /3(r) from optical observations alone. 
The filled circles in Figure 3 are the data of Fried & Illing- 
worth (1994). Both sets of observations lie entirely within 
the effective (half-light) radius, R,, = 104" = 8.57 kpc. 
For some reason the normalization is systematically dif- 
ferent in the two sets of optical observations. The Fried- 
Illingworth data have smaller error bars, less overall scat- 
ter and resemble in overall form the velocity dispersion 
profiles <t{R) of other well-observed ellipticals [e.g. Ger- 
hard et al. (1998) for NGC 6703]. The observational un- 
certainties apparent in the data sets shown in Figure 3 
make it difficult to determine with confidence any radial 
dependence of the stellar anisotropy. Consequently, we 
assume that the orbital anisotropy (3 in Equation (3) is 
constant over the range r > ro = 0.34 kpc. 

In Figure 3 we compare solutions of Equations (3) and 
(4) with the observed line of sight stellar velocity disper- 
sion in NGC 4472. Each solution curve a{R) is based 
on the linear gas temperature profile in Figure 2 and a 
particular value of /3. More information about the so- 
lutions shown in Figure 3 can be found in Table 1. Of 
the three solid lines designed to fit the Fried-IUingworth 
data, solution 2 with initial <Jr{ro) = 430 km s~^ and 
/3 = 0.7125 fits the data best. Solutions 1 and 3 show 
additional marginally acceptable fits although they fall re- 
spectively below and above the Fried-IUingworth data at 
smaller projected radii R. Evidently fi = 0.71 ±0.15 is the 
most likely anisotropy based on the gas temperature ob- 
servations. In some of our solutions of Equation (3) o",^(r) 
becomes unphysically negative at a radius r^ less than the 
outer limit of our integration in Equation (4), rt = 3Re] 
for these solutions we replace rt with r^. In solutions 4 - 
6, shown with dashed lines in Figure 3, we find that the 
anisotropy parameter must have values /3 = 0.63 ± 0.15 
to be in approximate agreement with the data of Bender 
et al. (1994). Once <Tr{ro) is selected for each of these 
dashed line solutions in Figure 3, the vertical position of 
the solution is extremely sensitive to /3, i.e., /3 must be 
specified to one part in 1000 to pass through the center of 
the Bender et al. data. One may doubt if Nature can be so 
discriminating. Nevertheless, for either set of optical data, 
the corresponding values of /3 are about twice as large as 
the f3{r) < 0.34 foimd by Kronawitter et al. (2000) from 
the Bender et al. data. 

2.2. Determining T{r) from /3(r) 

The preferred orbital anisotropy in NGC 4472 deter- 
mined by Kronawitter et al. (2000) from the Bender et al. 
observations can be accurately represented with 
PKir) = O.U9rlJ^^''^exp{-rkpc/2A72), 
which peaks at r w 1.65 kpc where f]^ ~ 0.34. In this sec- 
tion we assume that f3K{r) is correct and solve Equation 
(3) for c(r)^ and T{r). As before we consider only lin- 
ear temperature variations defined by two temperatures 
Ti and T2 evaluated respectively at ri = 2.18 kpc and 
r2 = 8.55 kpc. We seek values of Ti and T2 for which 
solutions a{R) of Equation (4) pass through or near the 
Bender et al. optical data. 



We begin with a set of solutions for which f3 = Pxir) 
and the temperature profile is constrained to agree with 
the observed linear T(r) in Figure 2. When £*(r) and 
p{r) are taken from Figure 1 as before, the only remaining 
undetermined parameter required to solve Equations (3) 
and (4) is (Tj.(^o)- The four short dashed lines in Figure 4 
show the line of sight velocity dispersion profile ^{R) for 
269.10 < (Jr{To) < 269.25 km s~^, corresponding to solu- 
tions 7 - 10 in Table 1. Since we are using the Kronawitter 
et al. /3k (r), it is appropriate to compare these solutions 
only with the data of Bender et al. from which (3k {r) was 
derived. Of these, the best-fitting solution is 9, but the 
overall fiat or positive slope [da{R)/dR > 0] at larger R 
does not resemble the a{R) profiles of most well-observed 
E galaxies. Since none of these solutions is a particularly 
good fit to the Bender et al. data, it appears that Pxir) 
found from optical data alone is somewhat inconsistent 
with the observed hot gas temperature profile. 

To explore alternate X-ray temperature profiles that 
may agree better with the cr(i?) of Bender et al., we con- 
sider solution 11 in which the hot gas is isothermal at the 
mean value of the two Buote temperatures in Figure 2, 
Ti = T2 = 1.137 X lO'^ K. This solution, shown with a dot- 
ted line in Figure 4, decreases slightly with R, in better 
agreement with data in i? < 2 kpc. Solutions with negative 
dT/dr may fit even better, but all known X-ray observa- 
tions for E galaxies have positive dT/dr for r < Re- Fi- 
nally, we consider solution 12 based on a temperature pro- 
file having the same slope as the observations in Figure 2, 
but with uniformly higher temperatures, Ti = 1.306 x 10^ 
K and T2 = 1.768 x lO'^ K. In solution 12 (solid line in 
Figure 4) the radially decreasing line of sight dispersion 
(T(i?) follows the Bender et al. data better than any other 
curve in Figure 4, possibly suggesting that the gas tem- 
perature should be higher than the observations in Figure 
2. Therefore, either the gas temperatures from X-ray ob- 
servations in Figure 2 are too low or PK{r) and the (t{R) 
data from which it was derived are inaccurate. 

3. SUMMARY AND FINAL REMARKS 

We have shown that Equations (3) and (4) establish a 
useful correspondence between the anisotropy of stellar or- 
bits in elliptical galaxies /3(r) and the temperature of the 
hot interstellar gas T{r). In practice, the /3 — T relation 
is based on the assumption of hydrostatic equilibrium for 
the hot gas and requires reliable observations of the stel- 
lar and gas densities and the stellar line of sight velocity 
dispersions. This connection between /3 and T is possible 
only if the galaxies are nearly spherical, as emphasized by 
Maggorian & Ballantyne (2002). 

If the gas temperature profile T(r) in an E galaxy is se- 
curely known from X-ray observations, P{r) can be deter- 
mined from the observed stellar velocity dispersion using 
Equations (3) and (4). In principle, this method of deter- 
mining /?(r) could be done at galactic radii much larger 
than Re where stellar observations of hi{R) may be diffi- 
cult but where observations of (j{R) and X-ray gas tem- 
peratures are generally reliable. Conversely, if [3{r) can 
be observed with high precision from optical observations, 
the hot gas temperature T{r) can be found from Equations 
(3) and (4) by using the gas density from X-ray observa- 
tions. This means of determining T(r) may be particularly 
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valuable in regions of the hot gas atmospheres where X- 
ray spectra indicate a two-tcmpcrature or limited multi- 
temperature thermal structure (e.g. Buote 2002b) that 
may complicate the determination of the average local gas 
temperature. 

As a representative example, we have used the (3- T re- 
lation to determine /3(r) in the bright elliptical galaxy 
NGC 4472 from the hot gas temperature profile T(r) 
found from X-ray observations. The resulting anisotropy, 
/? ^ 0.71 ± 0.15, indicates that the stellar orbits arc con- 
siderably more radially biased than implied by the smaller 
Pk{t) ^ 0.3 determined from optical observations alone 
(Kronawitter et al. 2000). These large /3 for NGC 4472 
are based on the Fried-IUingworth stellar velocity disper- 
sions a{R). Somewhat lower anisotropics. (3 = 0.63 ± 0.15 
are consistent with the velocity dispersion observations of 
Bender et al., but these /? still exceed the /3if(r) ^ 0.3 
found optically by Kronawitter et al. (2000). Alterna- 
tively, if the Kronawitter et al. /3 = /3_fc(T) and the Bender 
et al. (1994) data are correct, then the gas temperature in 
NGC 4472 may be higher than previously thought from X- 
ray observations and possibly have a different slope dT/dr. 
These conclusions apply only to the limited region r < Rg 



for which a{R) data currently exist for NGC 4472. 

It is clear, at least for NGC 4472, that currently avail- 
able observations are not sufficiently accurate to establish 
either f3{r) or T{r) with complete confidence. Accurate 
optical observations of (j{R) and stellar line profiles can 
be observed with 8-10 meter telescopes to several Re and 
T(r) is also often known at these large radii. In addition, 
it may be possible to determine gas temperature profiles 
with more accurate or prolonged X-ray observations than 
those in Figure 2. The /3(r) determined from X-ray emis- 
sion can serve to calibrate the accuracy of purely optical 
determinations or to provide a test for hydrostatic equilib- 
rium. To accomplish this, it would be desirable to obtain 
high quality stellar dispersion and gas temperature obser- 
vations of those galaxies such as NGC 4472, NGC 4649, 
NGC 1399, M87, etc. - that are luminous at both optical 
and X-ray frequencies. 



Studies of the evolution of hot gas in elliptical galaxies 
at UC Santa Cruz arc supported by NASA grants NAG 
5-8409 NAG 5-13275 and NSF grants AST-9802994 and 
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Fig. 1. — The observed and azimuthally averaged electron density n (cin^'"') profile in NGC 4472 as a function of radius normalized to the 
effective radius Re = 8.57 kpc (distance o! = 17 Mpc). The observations are from Einstein (Trinchieri, Fabbiano, & Canizares 1986) (filled 
circles) and ROSAT (Irwin & Sarazin 1996) (open circles). For the inner region we have Abel-inverted Chandra surface brightness data from 
Loewenstein et al. (2001) (open squares) and normalized them to jHCvious observations. The solid line is an analytic fit to the observations. 

The dashed line is the square root of the stellar luminosity density iV^(r) with arbitrary vertical normalization. 

Fig. 2. — X-ray observations of the gas temperature in NGC 4472 based on deprojected ROSAT observations of Buote (2000) {filled 
circles). Also shown with a filled square is the Chandra observation of the gas temperature at 0.082 kpc from Soldatenkov et al. (2003). Note 
that the central temperature is nonzero. The solid curve is a linear fit for T(r) that passes through all three observed temperatures. 

Fig. 3. — The line of sight stellar velocity dispersion as a function of projected radius. The observations are from Bender et al. (1994) (open 
circles) and Fried & lUingworth (1994) (filled circles). The lines show constant /? solutions based on the gas temperature T(r) variation shown 
in Figure 2. Each line type corresponds to several numbered solutions listed in Table 1 where values of <T(ro = 0.34 kpc) or cr(Re = 8.57 kpc) 
can be used to uniquely identify each line in the plot. 

Fig. 4. — The line of sight stellar velocity dispersion as a function of projected radius. The data arc the same as in Figure 3. The lines 
show solution curves based on the orbital anisotropy I3x(r) from Kronawitter et al. (2000). The numbers associated with each curve refer to 
solutions listed in Table 1. Each short dashed line can be uniquely identified with a solution in Table 1 by comparing values at a(Re) = 8.57 
kpc. 



Table 1 

SOLUTIONS OF EQUATIONS (3) AND (4) 

Solution Ti^ Gr{ro) a^(i?e) r/ u{ro) a{Re) 

(10^ K) (10^ K) (kms-i) (kms-i) (kpc) (km s'^) (km s'^) 



1 


0.6333 


.906 


1.368 


390. 


249. 


25.6 


282. 


172. 


2 


0.7125 


.906 


1.368 


430. 


253. 


23.7 


293. 


162. 


3 


0.7775 


.906 


1.368 


470. 


267. 


> 3i?e 


306. 


171. 


4 


0.4766 


.906 


1.368 


330. 


276. 


> 3i?e 


306. 


295. 


5 


0.63 


.906 


1.368 


390. 


290. 


> 3i?e 


318. 


270. 


6 


0.773 


.906 


1.368 


470. 


326. 


> 3Re 


334. 


245. 



13 



7 




.906 


1.368 


269.10 


205. 


> 3i?e 


247. 


210. 


8 


13k 


.906 


1.368 


269.15 


219. 


> 3i?e 


272. 


281. 


9 


Pk 


.906 


1.368 


269.20 


232. 


> 3Re 


295. 


337. 


10 


Pk 


.906 


1.368 


269.25 


244. 


> 3Re 


317. 


385. 


11 


Pk 


1.137 


1.137 


330.38 


254. 


> 3i?e 


315. 


303. 


12 


Pk 


1.306 


1.768 


332.9 


258. 


> 3Re 


316. 


295. 



^Pk refers to /3{r) from Kronawittcr ct al. (2000). 

''The linear temperature variation is defined by values Ti at ri = 2.18 and T2 at r2 = 8.55 kpc. 
'^The radial velocity dispersion (7r{r) goes to zero at r^. 
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